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Abstract 

We propose an improved algorithm for counting the number of Hamiltonian cycles in 
a directed graph. The basic idea of the method is sequential acceptance/rejection, which 
^ is successfully used in approximating the number of perfect matchings in dense bipartite 

graphs. As a consequence, a new ratio of the number of Hamiltonian cycles to the number 
^ ^ of 1-factors is proposed. Based on this ratio, we prove that our algorithm runs in expected 

[ time of 0{n^'^) for dense problems. This improves the Markov Chain Monte Carlo method, 

O^l ■ the most powerful existing method, a factor of at least ri*'^(logn)^ in running time. This 

class of dense problems is shown to be nontrivial in counting, in the sense that they are 
CN ■ #P-Complete. 
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1. Introduction 

A Hamiltonian cycle is a closed directed path that visits each vertex once and 
only once. In this paper we use digraph to denote directed graph. Counting the 
number of Hamiltonian cycles is a very challenging problem and has applications, 
for example, in quantum physics |4]. Many intractable counting problems have been 
added to the Valiant 's|20j list of ^^P-Complete, which is a natural correspondence of 
the concept NP-Complete for decision problems. Efficient approximating schemes 
called fully polynomial randomized approximation scheme(FPRAS) are naturally 
considered for the hard problems in counting. If M is the true value, a randomized 
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algorithm is called an FPRAS if it takes polynomial time of size of inputs, e and 
log (5^^ to obtain an output M. Here M is the approximation of M, satisfying 

P((l - e)M < M < (1 + e)M) >l-5. 

Due to the fact that the decision problem of whether a graph contains a Hamilto- 
nian cycle is NP-Complete, there would be no FPRAS for counting the Hamiltonian 
cycles for general graphs unless NP=RP. Thus the FPRAS for counting Hamilto- 
nian cycles are only possible for special or restricted graphs, for example, elementary 
recursive algorithms [T7j for random digraphs; Markov Chain Monte Carlo(MCMC) 
methods for dense undirected graphs [6j, for some random digraphs [8] and random 
regular graphs [7]. 

Sequential acceptance/rejection method is introduced by Huber [9j for counting 
the number of the perfect matchings in a dense regular bipartite graph. Recently 
the regularity requirement is removed [lOj. The primary tool used in the algorithm 
is the generalized Bregman's bound and the matrix scaling method. 

The MCMC algorithm presented for random digraphs in ^ can be naturally 
extended to dense digraphs. This algorithm is based on sampling 1-factors of the di- 
graphs and uses the self-reducing method [T3j to approximate the counting. Recently 
Bezakova et. al. present an algorithm that approximates the number of 1-factors in 
0(n''(log n)^) expected time, via an accelerating simulated annealing technique[l]. 

The ratio of the number of 1-factors to the number of Hamiltonian cycles is 
established to be 0(ni+i/(2"-3/2)) this paper provided that the digraph is an 
dense. Due to this ratio and Bezakova's results in [1], MCMC method[8j runs in an 
0(n^'^ (log n)'*) time when a > .85. Moreover, counting the number of Hamiltonian 
cycles in such digraphs is shown to be still #P-Complete. 

Our algorithm for counting Hamiltonian cycles is built on the acceptance/rejection 
algorithm in [10] while a different sequential sampling procedure is constructed to 
ensure that the approximating target is the number of Hamiltonian cycles. 

One of the remarkable advantages of acceptance/rejection method is that it sam- 
ples perfectly from a given set, which removes the sampling error when the MCMC 
method is adopted. Hence, our algorithm generates a weighted Hamiltonian cycle ex- 
actly according to its weight from the set of Hamiltonian cycles of a weighted digraph. 
In addition, this perfect sampling is only by-product when acceptance/rejection is 
used to approximate counting, which means the time used to sample a random 
Hamiltonian cycle can be used to approximate the number of the Hamiltonian cy- 
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cles without extra cost. The main result of this paper is summarized in the following. 

Theorem M. For any £, S E (0,1] and a G (.75,1], there exists a randomized 
approximation algorithm which provides an FPRAS for computing the number of 
Hamiltonian cycles of an dense digraphs. The same algorithm approximates the 
number of Hamiltonian cycles by a factor in [1 — e, 1 + e] with probability at least 
1 — S and has the complexity 0(n^-^+-^/(^"~^)+^/(^"~^-^)£~^ log{d~^)). In particular, 
when a > .85, the running time is bounded by 0{n^'^). 

The remainder of the paper is structured as follows. In Section 2 some basic 
definitions, notations and lemmas are presented. In Section 3 we describe the algo- 
rithm in details. Section 4 contributes to the complexity of the algorithm and the 
hardness of counting. Further discussion and conclusion are proposed in Section 5. 

2. Preliminaries 

Consider a simple weighted digraph G = {V, E) with the vertex set F = {1, ■ ■ ■ , n} 
and the edge set E. Each edge (i,j) G is endowed with a positive weight 
Wij. Let I S I denote the cardinality of any set S. The set of vertices point- 
ing to i is denoted by N~{i,G) = {j : {j,i) G E}, and similarly that out of i 
by N'^{i,G) = {j : {i,j) G E}. Indegrees and outdegrees of the vertex i are 
denoted by A-{i) = \N~{i,G)\ , A+{i) = \N+{i,G)\ respectively. Let A{i) = 
min(A~(i), A+(z)) and A = minjgy A(z). G is called an dense if A > an for an 
a > given. Let © denote the symmetric difference of two sets and [nj denote 
the maximum integer no more than n. A/B is used to denote the set by removing 
elements of B from the set A. With a little abuse of notation, / also denote the 
quotient of two numbers. A Hamiltonian cycle in G is represented by 

H = {ki,k2,--- ,k„,ki), 

where {ki, ^2, • • • , fc„} is a permutation of {1, • • • , n} such that {kn, ki) G E, and(/cj, 
/cj+i) e E, j = 1, - ■ ■ , n — 1. The length of a cycle or path is defined as the number 
of its edges that contains. 

An 1-factor is defined as a spanning directed subgraph of G in which indegrees 
and outdegrees of each vertex are all one. An example of an 1-factor is a spanning 
union of vertex disjoint directed cycles. Obviously, a Hamiltonian cycle is a special 
1-factor with only one cycle. The weight W{F) of an 1-factor F with edge set 
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{e G E}e^F is defined as W{F) = OeeF^e- The total weight W{S) of the set S 
of 1-factors are defined as W{S) = J^Fes^i^)- Let Wf{G) and Wh{G) denote 
the total weight of all the 1-factors and Hamiltonian cycles in G respectively. It is 
easy to see if Wij = 1 for all G E, then Wf{G) and Wh{G) are the number of 
1-factors and Hamiltonian cycles in G respectively. 

Let Ag be the adjacent matrix associated with G where AG{i,j) = Wij if 
G E and AG{i,j) = otherwise. For an n x n matrix A, where n is the 
order of A, we use notation Aij to denote the (n — 1) x (n — 1) matrix obtained 
from A by removing row i and column j. If there is no confusion, A^^j or (Aij) 
denotes the (n — l)x(n — 1) matrix obtained from A by, first permutating row i and 
row j and then removing row j and column j. Next we will define two quantities 
on the matrix Ag which are related to 1-factors and Hamiltonian cycles respectively. 

Definition 1. The permanent of an n x n matrix A = {A{i,j))nxn is 

n 

pev{A)=J2ll'^ih<^ii)), 

a i=l 

where a ranges over all the permutations of {1, • • • ,n}. 

Definition 2. The Hamilton of an n x n matrix A = (^(^,i))nxn is defined as 

ham(A) = A{ki, l)A{k2, h) ■ ■ ■ A{kn-i, kn-2)A{l, kn-i), 

{kl--- ,kn-l} 

where {ki • • • , kn-i} ranges over all the permutations of {2, • • • , n} when n >2, and 
ham(A) = A{1,1), ifn = 1. 

By the definition of permanent and Hamilton, it is not difficult to see that 
per(yl) > ham(yl) if the entries of A are all nonnegative. Suppose A = Ag- 
For any permutation (fci • • • of (2, • • • ,n), Aciljkn^i), Acikn-i, k„-2), 

AG{k2, kl), Aciki, 1) are the edge weight of the Hamiltonian cycle (1, kn-i, ■ ■ ■ ,ki,l) 
in G if and only if they are all positive. Therefore, we have 

Wh{G) = ham(^G)- 

Note that the diagonal entries of Ag arc all zero, and for any permutation a over 
{1, 2, • ■ ■ , n}, A{i, cr{i)) > 0, i = 1, • • • , n if and only if their corresponding edges in 
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G form an 1-f actor of G. Hence 

Wf{G) = per(AG). 

Next we present the Laplacian expansion formulas for the permanent and the 
Hamilton. 

Lemma 3. Let A = [A{i,j)) 

nxn be an nxn matrix. The permanent of empty 

matrix is set to be 1. Then 

n 

per(A) = ^A(i,l)per(Aii). 

i=l 

Lemma 4. Let A = {A{i, j))nxn be an n x n matrix, n > 2. Then 

n 

ham(^) = ^ A{i, 1) ham{A[i). 

i=2 

For the permanent, this expansion is well known. For the Hamilton, the formula 
is very similar and [17J proposes a combinatorial proof when each edge weight of 
the digraph is one. Regarding its importance in our algorithm, a proof in terms 
of matrix is presented below. We emphasize Lemma 4 is crucial in the sequential 
sampling procedure which is different from the one used in [10], and ensures our 
algorithm to approximate the number of Hamiltonian cycles. 

Proof of Lemma 4. We proceed to prove the lemma by induction on n, the 

order of the matrix. 

The case /c = 2 is trivial. 

Suppose Lemma 4 holds for k = n — 1. 

Consider k = n. Since 

n 

ham(A) = A{k2,i)A{k^,k2)--- A{kn-i,kn-2)A{l,kn-i), 

«=2 {fc2--- ,A:„_i} 

it is sufficient to show that 

ham(^-i) = ^ A{k2,i)- ■ ■ A{kn-i,kn-2)A{l,kn-i), 

{k2--- ,fc„-i} 
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for i = 2, • • • , n, where {k2 ■ ■ ■ , goes over all the permutations of {2 ■ ■ ■ ,n}/i. 

Considering the definition of A'-^, the row z — 1 of A'-^ is the first row of A except 
removing the first element, and 

A[^{k2-l,i-l) = A{k2,i), , 

A'-^{kn-i - 1, kn-2 - 1) = A{kn-i, kn-i), and 
- 1, K-i - 1) = ^(1, K-i)- 
By the hypothesis of the induction, the order of A^^ is n — 1, then 

ham(A-i) = ^ ^■i(/c2,i-l)----4^i(/e^-i,^n-2)^a(^-l>^n-i) 

= A;i(fc2-l,i-lKl(fcn-l-l,fen-2-l)^a(^-l>fcn-l-l) 
= ^ A{k2,i)---A{kn-l,kn-2)A{l,kn-l), 

{k2— ,kn-l} 

where {k2 • • • , and {k2 • • • , fcn-i} go over all the permutations of {1, • • • , n — 

l}/{i — 1} and {2, • • • ,n}/{i} respectively. This completes the proof of Lemma 4. 
□ 

Hamiltonian Recovery Let A = {A{i,j))nxn be an n x n positive matrix. The 
following procedure is applied to selecting elements from A (The first two steps are 
given explicitly). We call this procedure Selecting Hamiltonian Cycle(SHC for sim- 
plicity). 

Step 1. Let A^ = A. Choose a natural number I < ji < n, denote 7r(l) = ji 
and select ^^(7r(l), 1). 

Step 2. Let A^ = {Aj_^i)' . Choose a natural number 1 < j2 < n — 1, denote 
7r(2) = j2 and select A'^{Tr{2),l). 

Similarly A'', Tr{k) and A''{'K{k), l),l<k< n — 1 can be obtained in Step k itera- 
tively. Since A^ has only one entry, let A" = (^2^\i)' , 7r(n) = j„ = 1 and select 
A"(7r(n),l). 

By Lemma 4, the set of selected elements A^{Tr{k), 1), fc = 1, • • • , n, from the 
above procedure forms the edge weight of a Hamiltonian cycle in G if A = Aq- 
If 7r(l),7r(2), • • • ,7r(n) is given, we provide a simple algorithm to determine which 
Hamiltonian cycle in G is selected. This process is called Hamiltonian Recovery. 

The input of the algorithm is tt = (7r(l), 7r(2), • • • ,7r(n)). We illustrate how to 
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recover an entry in if 7r(l) is given. Let J^ii^j) be any entry in . Since 
= (^^(1)1 )' and recall the definition of which is obtained by, first 

permutating row 7r(l) and the first row and then removing the first row and first 
column. Hence, Hi = 7r(l) — 1 then (1, j + 1) is the position where A'^{i,j) lies of A^; 
otherwise {i + 1, j + 1) is the position where A'^{i,j) lies of A^. Hence, if the vector 
(7r(l), • • • , TT{k — 1)) is given from the SHC procedure, the position of A'^{7T{k), 1) in 
A^ can be found recursively by determining its position in A'^"^, then in A'^~^, and 
finally in A^. Since at each step of the SHC procedure an element is selected from 
the first column, A''{Tr{k), 1) must lie in column k of A^. 

If (1, /ci, • ■ ■ , kn-i, 1) is the corresponding Hamiltonian cycle of tt = (7r(l), 7r(2), 
••• ,7r(n)), then ki can be obtained from fcj+i since the element A{ki,ki+i) is se- 
lected in Step ki+i of the SHC procedure, or equivalently {ki, h^i) is the position of 
A'^'+i(7r(A;i+i), 1) in A, i = 1, 2, • ■ ■ ,n — 2. Obviously, kn-i = 7r(l). By this simple 
procedure, it takes 0{n?) time to recover all the positions of ^'^(7r(A;), 1), 1 < k <n. 
We present the recovery algorithm explicitly. 

Hamiltonian Recovery Algorithm 

Input : The vector (7r(l), 7r(2), • ■ ■ , 7r(n)). 
Output : A Hamiltonian cycle {l,ki, - ■ ■ , kn-i, 1). 
Step 1: Set kn-i = 7r(l); 

For z = n — 2 to 1 

Set a = 7r(A;j+i); 
For j = fcj+i to 2 

If a = Tr{j - 1) - 1; Set a = 1; 
Else Set a = a + 1; 
End; 

Set ki ^ cij 

End; 

Goto Step 2; 
Step 2: Output (1, fci, • ■ ■ , kn-i, 1). 

For simplicity, let HR^-k) denote the output of the Hamiltonian Recovery Algo- 
rithm when the input is tt = (7r(l), 7r(2), • ■ • , 7r(n)). 

3. Algorithms for Counting 
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One main tool in our algorithm is a generalized version of Bregman's bound for 
the permanent below, which generalized an inequality of Soul [19] and proved in 
[lOj . For more application of other generalization of Bregman's bound for designing 
new algorithms or improving efficiency of algorithms, we refer to [14: \ I18 j . Let 

( r + (1/2) logr + e - 1, r>l 

9ir) = { 

il + (e-l)r, re [0,1]. 

Lemma 5. ([lOj) Let A be an n x n matrix with entries in [0,1]. Denote r[i) 

n 

the sum of row i of A. Define Br{A) = Y\{g{r{i))/e), then 

i=l 

n 

Br{A)>Y,MiA)Br{Ai). 

i=l 

In particular, by Lemma 3, per(j4) < Br (A). 

Chernoff's bound is useful in our algorithm, and one form of that is given 
bellow [E]. 

Lemma 6. Let xi,X2,--- ,xt be identical independent distributed(i.i.d.) Bernoulli 
random variables with P{xi = 1) = p and P{xi =0) = 1 — p, p > then for any 
< e < 2e - 1, 

t 

P{\^Xi-tp\> etp) <e-'P'^/^. 

i=l 

For simplicity, in this section we only consider the digraph G with all edge weight 
equalling one. Hence the adjacent matrix Aq is a 0-1 matrix and ham(A(;) is the 
number of Hamiltonian cycles in G. G is also restricted to be an dense, a > .75. 
It is known [2] that if G is .5n dense, G must contain a Hamiltonian cycle and the 
proof can be easily modified to give an O(n^) algorithm to construct a Hamiltonian 
cycle. Hence ham(AG) > 1. By the definition of Hamilton, if we change any zeros 
in Ag to 7 = (e/3)((n — l)!)^"*^, ham(^G') increases by at most a factor of 1 + e/3. 

Now we introduce the basic idea of acceptance/rejection method for the counting 
problem. Suppose S* is a large set and each element in it with positive weight. The 
target is to approximate the total weight of all the elements in S. First select a 
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suitable large M such that M > YlteS ^(^)- The main idea of acceptance/rejection 
method for approximation is to design a procedure to sample a random element x 
from the set S with the successful probability P{x = a) = where w{a) is the 
weight of a € S, and failing probability P{x ^ S) = 1 — ^^^^§f-^- At each time, if a 
random element a is successfully selected from S, we say acceptance or a is accepted, 
and if no element is selected from S, we say rejection. Hence, at each time the 

probability of acceptance is — and probability of rejection 1 . With 

some fundamental statistical knowledge, the total weight of S can be approximated 
by multipling M and the ratio of acceptance over all the samplings. For our purpose, 
generalized Bregman's bound in Lemma 5 provides such a suitable large M, and self- 
reducing method for counting Hamiltonian cycles naturally proposes such a sampling 
procedure, which is sequential sampling procedure guaranteed by Lemma 4. For 
more details about sequential acceptance/rejection method, we refer to [lO] • 

In order to make use of the generalized Bregman's bound in Lemma 5, before 
resuming the acceptance/rejection algorithm, we need to scale the matrix Aq to 
nearly be doubly stochastic and make each entry in [0.1]|15j. Hence the algorithm 
has two phases. 

Sub Algorithm I. Scale Matrix 

Input : Ac, e 
Output : X, Y, Z, C 

Step 1: Set AdiJ) = (e/3)((n - 1)1)-^ if AdiJ) = for all i, j, goto Step 2; 

Step 2: Using matrix scaling to find diagonal matrix X, Y such that the row and 

column sums of B = XAqY in (1 — 1 + (.l)n~^), goto Step 3; 

Step 3: Let Z he a diagonal matrix with Z{i,i) = min^ B{i,j)~^ for i = 1, • ■ ■ n, 

goto Step 4; 

Step 4: C = ZB. 

After matrix scaling, matrix C satisfies the requirement of generalized Bregman's 
bound. Sequential acceptance/rejection method can be used to estimate ham(C). 
Note that the matrix C is corresponding to a weighted digraph denoted by Gc- 

Sub Algorithm II. Approximating Hamilton via Acceptance/Rejection 

Input : X, Y, Z, C, e, 6 N. 

Output : TCi, ■ ■ ■ jTi-s] ham(ylG) the estimator of ham(^G). 
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step 5: Set t = m{e/2)-^ log{6-^), I = Y[{X{i,i)Y{i,i)Z{i,i)), D = C, k = 

1=1 

and s = 0, goto Step 6; 
Step 6: Set r = order of D; 
If r = 1; 

Set p(l) = D/Br{D) and p(0) = 1 - p(l); 

Choose / from {0, 1} according to P{I = i) = p{i), i = 0, 1; 

If / > 0; Set 7r(n) = I, s = s + 1, k = k + 1 and Hs = HR{tt); 

If k < t: Set D = C, goto Step 6; Otherwise goto Step 7; 
Else / = 0; Set A; = A; + 1; 

If A; < i; Set D = C, goto Step 6; Otherwise goto Step 7; 

Else r > 1; 

Set p(i) = D{i,l)Br{D'-^)/Br{D) for i = 2, • • • ,r andp(O) = l-Ep(i); 
Choose / from {0, 2, 3, • • • , r} according to P{I = i) = i = 0, 2, • • • , r; 
If / > 0; Set 7r(n + 1 - r) = I and D = D'j^, goto Step 6; 
Else / = 0; Set /c = /c + 1 

If A; < t; Set D = C, goto Step 6; Otherwise goto Step 7; 
Step 7: ham(^G) = l-^st-^Br{C). 

The procedure of samphng elements in Step 6 is the same as SHC proccdme 
except selecting an element with certain probability or rejection when I = is se- 
lected. The output Hi, 1 < i < s, is accepted by the algorithm. 

Theorem 7. Let Tti,--- ,7is CL^d ham(^G) be the output of Sub Algorithm II. 
If we set N = Br{C) / ham{C) in the same algorithm, and let H be a random vari- 
able recovered from a random vr of Sub Algorithms II and S denote the set of all the 
possible accepted hamiltonian cycles, then 

p{H = ni\Hes) = w(ni)/WH{Gc) 

and 

P((l - e) ham{AG) < hamiAo) <(! + £) ham(AG)) >l-d. 

Proof. First, we check p(0) > at each level of Step 6, which guarantees the 
proceeding of the algorithm. By the definition of D'-^ and Da, obviously, Br{D'--^) = 
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Br{Dii). Using Lemma 5, it is easy to see 

n n 

D{i, l)Sr(Ai) = J2 l)Sr(Ai) 

i=2 i=2 
n 

<J2D{i,l)Br{Da) 

i=l 

< Br{D). 

Hence p(0) > 0. Suppose Hi = HR{j), j = (ji,--- ,jn)- Following the path in 
which Til is selected, and using the notation in SHC procedure, then C*"^^ = (Cj.^^)', 
i = 1, - ■ ■ ,n — 1 and = C, we have 



P{7rik)=jk) = 



C>^{jk,l)Br{{Cl,)) 
Br{C^) 



where A; = 1, 2, • • • , n - 1 and P(7r(n) = j„) = ^^^^.^'j ■ 

Since the selection at each level in Step 6 is independent of the other, the prob- 
ability of selecting Hi is the telescoping product. Noting that = (Cj.i)', 
i = 1, - ■ ■ ,n — 1 and = C, then 

Since each Hamiltonian cycle in Gc can be accepted with certain probability 
proportional to its weight, the acceptance set S is the set of all the Hamiltonian 
cycles in Gc- Then 

Wh{Gc) 

r{n = n) = 

neGc 

Hence, 

In Sub Algorithm H, let x^, 1 < A; < i, denote the indication function of ac- 
ceptance or rejection in Step 6, that is, Xfe = 1 if a Hamiltonian cycle is accepted 
and Xk = otherwise. Obviously, x^, ^ < k < t, are i.i.d. Bernoulli random 
variables with P{xi = 1) = p = WH{Gc)/Br{C) = ham{C)/Br{C). Let be 
the matrix obtained in Step 1 of Sub Algorithm I. Hence, by Lemma 6 and noting 
t = AN{s/2)~^ log{S~^), where N = Br(C)/ ham(C), a simple calculation shows 

P((l - e/2) ham(A^) < hain(^G) < (1 + e/2) ham(A^)) > 1 - 6. 



W{Hi) 
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Noting ham(Ac) < ham(j4g) < (1 + e/3) ham^Ac), thus the proof completes. □ 



4. Complexity and Hardness of Counting 

4.1. Complexity of the Algorithm 

Due to ellipsoid method|15j. the running time of matrix scaling is 0(n^ log n). 
So the complexity of Sub Algorithm I is 0(n^ logn). 

The time of repeating Step 6 in Sub Algorithm II is t = 0{Br{C) / ham(C)), and 
for each time the running time is O(n^), hence, the complexity of Sub Algorithm 
II is 0{n^ * t) = 0(n^i?r(C)/ ham(C)), where e~'^log5~^ has been put into the 
term O(-) for simplicity. As we know, the Hamiltonian Recovery Algorithm takes 
O(n^) time. After removing Hamiltonian Recovery procedure, the total running 
time of Sub Algorithm II is still 0{n^ * t), thus if approximating ham(AG') is only 
the purpose, outputting the Hamiltonian cycle is the byproduct of Sub Algorithm 
II. 

If the digraph G is an dense, a > .5, an important result given by Huber[10] is 

Br{C)/per{C) = ©(n-'^+'^/^^a-i))^ 

Note that 



per(C) i=i per(^; 



ham(C) fl^x{t,i)Y{i,i)Z{i,i))ham{Ag ham(A^) 

i=l 

If the digraph G is at least .5n dense, then changing any zeros in Aq to en~^ increases 
pei{AG) by at most a factor of 1 + e [11]. Then 

per(C) _ per(^^) ^ (1 + g/3) per(^G) ^ ^^_per(^G) , 



ham(C) ham(j4^) ham(AG') ham(j4G 

Hence, the total running time of our algorithm is 

^,4, 2 Br{C) , ^,4, 2Br{C) per(C) , 



ham(C)' ' ^ per(C)ham(C) 
0{n'^ logn + n 



4i„„„ , ^1.5+.5/(2a-l) Per(^G) 



(1) 



ham(^G) 

Now we present combinatorial argument on the bound of j^^^^^j (Recall Aq is a 
0-1 matrix and all the edge weight of G equals one). The methodology is analogous 
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to the approach for undirected graphs given by Dyer et.aL[B]. 



Lemma 8. ([6]) Let n be a natural number and (5 a positive number. Let ko 
max([/31ognJ, 1) and g{k) = n^k\{[3\ogn)~^ , define 



f{k) 



' g{k), k<ko 
_ g{ko), k > ko. 
Then f{k — 1) > {/3 log n)k^^ f{k); and f{k) > 1 for any k. 



Proof. Uk< ko, f{k - 1) = g{k - 1) = {p log n)k-'^g{k) = {piogn)k-'^ f{k); 
If k > ko, then piogn/k < 1. Hence 

f{k - 1) = g{ko) > {(3logn)k-^g{ko) = {(3 log n)k-^ f{k). 

Thus /(/c) > f{ko), we have 



Theorem 9. Suppose a £ (.75,1]. Let G be an an dense digraph and the 
set of 1 -factors in G containing exactly k cycles, 1 < k < [n/2\. Note that Fi is the 
set of Hamiltonian cycles in G. Let F = |J^ Fk- Then 

\F\ 1 + 1/(2q:-1.5)n 



1^1 1 

With this theorem, we prove the main result of this paper Theorem M. 

Proof of Theorem M. By theorem 9, since = per(AG)/ ham(^G'), and 

noting (1), therefore Theorem M fohows immediately. 

Now we proceed to prove Theorem 9. 

Proof of Theorem 9. We construct a new weighted digraph ^ = {F,K). K 
is defined as follows. 

K = {(E,e'):E £Fk,E' eF^,,k <k and E®e'^'CI}, 
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Figure 1: The symmetric differences Ci = E ® E 




Figure 2: E is obtained by coalescing two cycles of E into a single cycle 

where C4 is a graph with four vertices and four edges, in which two vertices have 
indegrees two, outdegrees zero, and the other two vertices have indegrees zero, out- 
degrees two(See Figure 1). The four edges belong to E and E' alternatively. To 
avoid the confusion with vertices and edges in G, we call the nodes and arcs cor- 
responding to F and K in ^. Observe also that if {E,E') G if is an arc of 
can be obtained from E by deleting two edges and adding two others, and that this 
operation can decrease the number of cycles by one(See Figure 2). Hence every arc 
{E, E') is directed from a node E in some Fj^ to a node E' in -Fj._i. 

The proof strategy is to define a positive weight function w on the arcs set K 
such that the total weight of arcs leaving each node E e F/Fi is at least one greater 
than the total weight of arcs entering E. Denote w'^{S) and w'~{S) the total weight 
leaving and entering a node set 5 in respectively, the strategy ensures 

w-{Fk) + \Fk\ = Yl ("'"(^) + 1) < E ^^(^) = ^^i^k), k>2. 

E&F^ E&Fk 

Hence, 

w-{Fi) = w+{F2) = - w-{Fk)) > ^ \Fk\ = \F/Fi\. 

k>2 k>2 
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Let g = maxE&Fi w (E). Since w (Fi) = J2e£Fi ^ i^) — then 

\F\/\F,\<g+l. 

The weight function w : K ^ defined as follows. For any arc {E, E ) with 
E' G Ffe, we know E' is obtained by coalescing two cycles of E, and suppose the 
length of these two cycles are li and /2, then define w{E,E') = (l^^ + l2^)f{k), 
where f{k) is defined as in Lemma 8. Then we have the following two claims. 

Claim 1. For any E eFk, k>2, w+{E) > (2a - 1.5)n/3/(fc) logn + 2. 

Claim 2. For any E € F^, k >1, w~{E) < nlognf{k). 

By these two claims, set /? = 1/(2q;— 1.5). Then for E € F^, k >2,we have w~^{E)~ 
w-{E) > 2 > 1 and 5 = m&yLE&FiW^ {E) < nlogn/(l) < (2a - 1.5)n^+i/(2"-^-^). 
Hence < g + l = 0{v}'^^/^'^°'~^-^^), which completes the proof. □ 

Proof of Claim 1. Let E € Fj. be an 1-factor with k cycles 71, • • • ,7^., of lengths 
ni,n2, ■ ■ ■ k > 2. We proceed to bound w~^{E). To show the lower bound of 
w~^{E), we need to count the number of arcs leaving E. Suppose {E, ) to be such 
an arc. Let 7 = E®E' , E' G -Ffe-i, be the form (x, x' ,y, y'), where (x, x'), {y, y') & E 
and {y,x'), {x,y') € E' . 

First, we estimate the number of C4-type cycles 7 for which {x,x ) is contained 
in a particular cycle 7j G E. We say that 7 is rooted at 7.;. Assume, for a moment, 
that the vertices x.x is chosen. There are at least an — {ui — 1) ways to extend 
the path first to y then to y since the indegrees of x' is at least an. Denote Y' the 
set of all vertices y reachable. Recall N^{x,G) is the set of neighbors x points to. 
Thus the number of ways of completing a C4-type cycle {x, x' ,y, y) is at least 

\N+{x,G)\ + \y'\ - n<an+{an - {rn -l))-n 
= 2an — rii — n + 1. 

Hence the total number of C4-type cycles rooted at 7^ is at least ni{2an — rii — n + 1) . 

We are now poised to bound w~^{E). Each arc {E,E') defined by a C4-type 7 
rooted at 7^ has weight at least n^^f{k — 1), which, by Lemma 8, bounded below 
by {piogn){kni)-'f{k), Thus 
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E':{E,E')eK 

k 

— ^ ni{2an — rii — n + l)n^^f{k — 1) 

i=l 
k 

= ^(2an -Hi- n)f{k - 1) + kf{k - 1) 
j=i 

> {2akn - n - kn) {(3 log n) k'^ f{k) + kf{k - 1) 
= (2a-l/A;-l)(^logn)/(A;)n + A;/(A;-l) 

> (2a- 1.5)(/31ogn)/(A;)n + 2. 

For the first inequality, it seems we have overcounted the weight, we explain 
the reason. When {x,x') is rooted at -fi and {y,y') lies in some 'jj if we extends 
{x,x ) to {y,y ) to complete a C4-type cycle 7 = {x,x' ,y,y'), the contribution to 
the weight is only n~^f{k — 1) in the above inequality. Similarly, when (x,x') is 
rooted at the same position as {y,y') in 7^ and {y,y') lies in the same position as 
{x,x') in 7j, the contribution to the weight is nj^f{k — 1). Plus these two weight, 
{nj^ + n^^)f{k — 1) is exactly w{E,E') needed to be considered by the definition 
of w, where E ® E' = ^. Hence, though each C4 cycle is counted twice, the weight 
not. The last inequality follows immediately from k>2 and f{k — 1) > 1. □ 

Proof of Claim 2. For each E G F^, we now proceed to bound w^{E). Let 
{E',E) be an arc in K. It is straightforward to verify that the C4-typc 7 = 
{x,x' ,y,y') = E ® E' must contain two edges and {y,y') from a single 7^ 

of E, and {y,x'),{x,y') G E' . Removing these two edges from 7j leaves a double of 
simple paths of lengths p — 1 and g — 1, where p,q > 2. For the case p ^ q there 
are at most rii ways such that 7^ © 7 is a pair of cycles with length p and q, and 
ni/2 ways such that 7^ © 7 is a pair of cycles with length p and q for the case p = q. 
Noting both cases happen when 7j is contained in a complete sub digraph of G or G 
is a complete digraph (Complete digraph is defined as such a digraph that any two 
distinct vertices have edges pointing to each other). Hence 
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w~{E)= w{e',E) 

E':{E',E)eK 




p+q=ni p-\-q=ni 



p+q=ni 




k ni~2 

1=1 p=2 ^ 

k 

< (k) log Ui 

i=l 

< nlognf{k). □ 



4.2. Hardness of Counting Hamiltonian cycles in dense digraphs 

We first declare the notation related to undirected graphs only appears in this 
subsection and the notation related to digraph is the same as that in the previ- 
ous sections. Our reduction comes from the undirected graph, hence notation for 
undirected graphs is needed. Let G be a simple undirected graph with vertices 
{1,2, ■ ■ ■ ,n}, where n > 3. The definition of a Hamiltonian cycle of an undirected 
graph is a closed undirected path that visits each vertex once and only once. We use 
the notation mim-2, ■ ■ ■ ninmi to denote a Hamiltonian cycle in an undirected graph 
(recall (mi,m2, ■ • • ,m„,mi) denotes a Hamiltonian cycle in digraphs). The degree 
of a vertex in an undirected graph G is defined as the number of its neighbors. Let 
#HC and #DHC be the problem of counting the number of Hamiltonian cycles in 
undirected and directed graphs respectively. Now define a symmetric digraph G' 
corresponding to an undirected graph G by replacing each edge of G with two 
directed edges and Let Hq and Hq' denote the set of the Hamiltonian 

cycles in G and G' respectively. Let V{Ho') denote the power set of Hq'- We will 
prove the number of Hamiltonian cycles in an undirected graphs equals half of the 
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number of Hamiltonian cycles in its corresponding symmetric digraphs. 

Lemma 10. ([6]) i^HC is jj^P- Complete, even when restricted to graphs G of min- 
imun degree at least (1 — e)n, for any e > 

Lemma 11. Let Ti = mi • • • m„m-i he a Hamiltonian cycle in Hq- Then there 
are at least two Hamiltonian cycles {mi, ■ ■ ■ ,m„,mi) and {m\,mn, ■ ■ ■ ,Tni) in He ■ 
Define a map f from Hq to V{Hqi) as follows: 

ip{TC) = {(mi, • • • ,m„,mi), (mi,mn, • • • ,mi)}. 

Let Irmp denote the image set of the map ip, and let 7i' = m'^ ■ ■ ■ m'^m[ be a different 
Hamiltonian cycle from TL in Hq . Then 

ip{n) n if{n') = and Ulmy? = Hg'. 

Proof. Due to the symmetry of the digraph G\ and noting n > 3, for any Hamil- 
tonian cycle (mi,-- - ,m„,mi) in Hq', there must be a different Hamiltonian cy- 
cle (mi,m„,-- - ,mi) in Hqi. These two Hamiltonian cycles obviously has a pre- 
imagine, the Hamiltonian cycle mi - - - m^mi in Hq. Note (mi,-- - ,m„,mi) is in 
(/^(mi - - - m^mi). Hence, Ulrmp D Hql Obviously, iMrmp C Hql Therefore 

Ulnup = Hqi. 

Suppose there are two different Hamiltonian cycles Ti = mi ■ ■ ■ m„mi and Ti' = 
m'l ■ ■ -m'^m'i in Hq. Let N-}-i{mi) denote two neighbor vertices of vertex mi in Ti. 
7i and 7i' are different if and only if there exits a vertex {mj}={m^} such that 
N'H.i'^i) 7^ -^H'(^j)- Hence (mi,-- - ,m„,mi) is different from (m'^,-- - ,m^,m'^) 
and (m']^, m^, • • • , m2, m'^), that is (mi, - - - , m„, mi) ^ fCH'). Similarly, {mi,mn, - - - , mi) 
^ ip{n'). Hence ip{n) n ip{n') =0. □ 

Theorem 12. f^DHC is jfP-Complete, even when the digraph is (1 — 7)n dense, 
< 7 < .5. 

Proof. Lemma 11 shows the number of Hamiltonian cycles in an undirected graph 
is half of the number of Hamiltonian cycles in its corresponding symmetric digraph. 
Hence by Lemma 10, #DHC in (1 — 7)77, dense digraphs is #P-Complete, for any 
< 7 < .5. □ 
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5. Conclusions and Discussions 

The results in this paper show that for relatively dense digraphs, approximating 
the number of Hamiltonian cycles or generating weighted Hamiltonian cycles exactly 
from their correct distribution can be accomplished in (9(7^2. 5+.5/(2«-i)+2/{4a-3)-j 
time. This is an improvement in running time by a factor of n^'^(logn)^ for .85n 
dense digraphs. Counting the number of Hamiltonian cycles in such digraphs is 
shown to be #P-Complete. 

Estimating the Hamilton of a 0-1 matrix to within a factor of 1 + e with proba- 
bility at least 1 — S, the running time is 

Q(^2.5+.5/(2a-l)+l/(2a-1.5)^-2i^g(^-l))_ 

It is known [2] that 0.5n dense digraphs contain Hamiltonian cycles. Our algo- 
rithm presented in this paper is shown to be an FPRAS for 0.75n dense problems. 
Hence a gap still remains. We can extend the definition C4 in the proof of Theorem 
9, as shown by Figure 1. Similarly that can also be done to Cg, Cg. However it 
seems unlikely to obtain any better bounds than that by C4 in this way. This gap 
is left open here. 
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